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We present an L-A pair for the Apel'rot case of a heavy rigid 3-dimensional body. Using it 
we give an algebro-geometric integration procedure. Generalizing this L-A pair, we obtain a new 
completely intcgrablc case of the Euler-Poisson equations in dimension four. Explicit formulae for 
integrals which are in involution are given. This system is a counterexample to one well known 
Ratiu's theorem. Corrected version of this classification theorem is proved. 

1. Introduction 

The rotations of a heavy rigid body fixed at a point are described by the Euler- 
Poisson equations. It is well known that these equations are integrable in Euler, 
Lagrange and Kovalevskaya cases. For a fixed value of one of the integrals, there 
are additional integrable cases, for example Goryachev-Chaplygin and Apel'rot (see 
[1], [5]). L-A pairs were known for all these problems except the last one (see [4], 
[7], [10], [11]). In this paper, in section 2, we present an L-A pair for the Apel'rot 
case. (Almost the same L-A pair serves in the Apel'rot gyrostat case.) Using 
it, we give an algebro-geometric integration procedure in section 3. The spectral 
curve is reducible and it consists of a sphere and a torus. The evolution of the 
pole divisor of the eigen-function of the L matrix is linearized on the Jacobian of 
the torus. It leeds to the solution of the initial problem up to rotation. For a 
complete solution one needs one step more: to integrate the Riccati equation with 
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double-periodical coefficients. Generalizing this L-A pair, we obtain, in section 4, 
a new completely integrable system of a four-dimensional rigid body motion. This 
system is integrable even without reduction to any invariant submanifold. Explicit 
formulae for integrals wich are in involution are given. 

It turns out that the L-A pair, which we use, is of the form analysed by Ratiu and 
van Moerbeke in the study of the generalized Lagrange case and of the completely 
symmetrical case in [9], [10]. Our system is a counterexample to Ratiu's well known 
theorem (see [9], [13]). This theorem claims that the Lagrange and the completely 
symmetric case have exhausted the list of the Euler-Poisson equations which are 
equivalent to an L-A pair of that form. We conclude the paper by giving a corrected 
version of this classification theorem. 



2. Apel'rot case 

The equations of rotation of a heavy rigid body fixed at a point in a moving 
frame are: 

^M=Mxw + 7xf c , ^7 = 7X(i, (1) 

where uj is the vector of the angular velocity, M — Iuj is the kinetic momentum, 
I is the inertia operator, 7 is the unit vertical vector, and fc — {xo,Uo, z o) is the 
radius vector of the mass center according to the fixed point. One can choose the 
moving frame such that the inertia operator is diagonal, / = diag(Ii, I2, ^3)- The 
equations (1) have the following three integrals of motion (see [3]): 

^1 = 7j(M<^) + (l,rc), 
*2 = <7,7>(= 1), 
F 3 = (Jw,7). 

For complete integrability, we need one integral more. In 1894 Apel'rot noticed 
(see [1], [5]) that under the additional conditions: 

*) 2/o = 

(2) 

it) xoy/h(h ~ h) + zoVWi -h) = 
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the equations (1) became integrable on the hypersurface given with: 

hxouji + I 3 z q uj 3 = 0. (3) 

A geometric interpretation of these conditions was given by Zhukovski [6], [14]. 

Using standard isomorphism between Lie algebras R 3 and so(3) (notation: A e 
R 3 A e so(3)), (see [2], [10]), the equations (1) can be transformed to the 
following equations on so(3) x so(3): 

^ = [MM + h,r c ], ^ = hM (4) 

Theorem 1. Under the Apel'rot conditions (2), the equations (4) are equivalent 
to: 

^(A 2 C + AM + 7 ) = [\ 2 C + \M + -f,uj + \r c }, (5) 
where X is a spectral parameter and C € so(3) is a constant matrix given by 

C = hr c - (6) 

Proof. The equations (4) are equivalent to (5) if and only if 

[C,r c ]=0, [C,Lo} + [M,r c ] = 0. 

The first equation is satisfied by (6). The second one can be reduced to the equation: 

(h - I 2 )uiz + (h - h)u 3 x a = 0. (7) 

But (7) is a consequence of the Apel'rot conditions (2) and the equation of the 
hypersurface (3). 



Note. A gyrostat is a system which consists of two rigid bodies. The second one, 
the gyroscope, is axialy-symmetric and it rotates with a constant angular velocity 
about its symmetry axis which is fixed in the first rigid body. The equations of 
motion of that system are of the form (4) with 

M = Ilo + P 
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where P is a constant skew-symmetric matrix, which represents the kinetic mo- 
mentum of the gyroscope. Sretenskiy in [12] proved that, if Apel'rot conditions 
(2) are satisfied, and P 2 = 0, then the equations of motion are integrable on the 
hypersurface given by: 

(h - h)uiz Q + (I 2 - h)usXo - (P 3 x - P\Zq) = 0. 

It can easely be proved that: 

L{\) = \ 2 C + \M + -(■ A(X)=oj + Xr c , 

form the L-A pair for the Apel'rot gyrostat. 

An L-A pair of the form 

L(X) = X 2 C + XM + 7, A(\)=iu + \r c (8) 

was used in [10] for the Lagrange case (of course, the matrices rc and C from [10] 
are different from the matrices rc and C in this paper. Also, the conditions on the 
inertia operator are not the same). 

3. Integration of the Apel'rot case 

The next steps in algebro-geometric integration of the Apel'rot case follow the 
procedure from [10]. 

As usual, in a direct problem one corresponds to L(A) some algebro-geometric 
data with simple evolution in time. After the integration of this evolution, in the 
inverse part, one has to reconstruct the matrix elements of L(X), starting from the 
algebro-geometric data as known functions of time. The final part is the integration 
of the equations of motion of the initial problem by use of the known matrix entries 
of L(X). A specific characteristic of the problem presented here, is that the algebro- 
geometric procedure gives a solution of the equations of motion of the Apel'rot case 
only up to rotations. For the full solution of the problem, we need one integration 
more to determine the angle variable (f> x = argx - see theorem 5 bellow. 
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Let: 



The matrix 



a 



x 



V4 + 4 



(9) 



U = 



a 



i§_ _§_-[ 
— 

P V2 V2 



diagonalises matrix C. In the new basis, the matrix L(X) = CA 2 + MX + 7 trans- 
formes into: 

L = U^LU = 



A iA* 
-A* -Q, 

iA n 



where: 



A = y + Xx; A* = y + Xx; 
y = -^=(/?7i - 0:73 - 172); x = -^=(pMi ~ aM 3 - iM 2 ) 
Sl = -i [aidX 2 + Mi A + 71) + /5(^ 3 A 2 + M 3 X + 73)] 
= -i [aidX 2 + 7l ) + (3(C 3 X 2 + 73)] 
and C\ = I 2 x ; C 3 = I 2 z , using (3). 
The spectral curve is defined by: 



(10) 



where 



T : p(n, X) := det(L(A) - fj,E) = 0, 



p{ji, X) = -n(p 2 - fl 2 + 2AA*). 



(11) 



The spectral curve T is reducible and it consists of two components: the sphere 
Fi given by /x = 0, and the torus T 2 : 



H 2 = n 2 - 2AA* =: P 4 (X). 



(12) 



The coefficients of the spectral polynomial (12) are integrals of the motion. Sub- 
stituting expressions (10) in (11), the following form of the equation of the spectral 
curve is obtained: 



p(ji, A) = -ii(/x 2 + AX A + BX 3 + DX 2 +EX + F), 



(13) 
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where 

A = I 2 2 (xl + zl) 

B = 2I 2 (M lXo + M 3 z )(= 0) 

/ M? Ml Ml \ 
D = 2I 2 ^ + ^ + ^ +X0 , 1 + Z0 ^ (14) 

E = 2(Mi7i + Af 2 72 + M373) 

F = ll+ll+ll{=l) 
Using the Apel'rot conditions, the expression for D can be simplified: 

„ „ r (Ml; Ml Ml 



, 2/i 2/ 2 2/3 

This is the energy integral. So, in this case the L-A pair (5) gives four integrals of 
the motion, which is enough for integrability. 

Let (/1, f 2 , fs) T denote an eigenvector of the matrix L(X), which corresponds to 
the eigenvalue fi. Fix a normalizing condition /1 = 1. Then: 

A* iA 

H = -7T— , h = -o • (15) 

The restrictions to T 2 are: 

_ O + it 

h ~ 2A~ ' h ~ 2iA* ' 
The relation between the divisors on T 2 : 

(/ 2 ) + (/a) = 0. 

is a consequence of / 2 • / 3 = -|. 

From (12) the divisors of A and it on r 2 are: 

(A) = -Pi -P2 + R1+R2 

(it) - -2Pj - 2P 2 + ( M ) , deg(it) = 4 

where Pi and P 2 are the points on T 2 over A = 00. From (12), we obtain the 
developments in neighborhoods of Pi and P 2 : 



il 2 \/x\ + Zq + 0(X 2 ), in a neighborhood of Pi 



X 2 1 —U 2 ^Xq + Zq + 0(A 2 ), in a neighborhood of P 2 . 
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From (15), it follows that f 2 has a simple pole in P\ and a simple zero in P 2 ; 
has a simple zero in Pi and a simple pole in P 2 . In the affine part, f 2 has a simple 
pole in v given by 

y 

v : A = 0, fl + ri = 0, i.e. v\ ~ , v u = —ft U=-i • 

x 

It has a simple zero in v = (v\, D^), where v\ denotes the complex conjugate of v\. 
The previous consideration gives us: 

Lemma 1. The divisors of f 2 and / 3 on T 2 are: 

(/a) = Pi - Pi + v - v 
Now we are going to analyse the converse problem. Suppose the evolution in 
time of the point v is known (see theorem 4 bellow). For reconstructing the matrix 
L(X), one needs x = |x|e larga: , y = \y\e iazsy . 

Theorem 2. The point i/£f 2 and the initial conditions for M and 7 determine 
\x\, \y\ and argy — argx, where x and y are given by (10). 

Proof. Let v have coordinates and v\. From the expressions for v 

y o 1 

x 

and expression (10) for ft, it follows: 

aji + /?73 = -ivy, - hv\\[x\ + ~z$- (16) 

From (10) and (16) we can determine as a polynomial of A: 

= -i[h\ 2 Jx~l + zl- ivy - hv\sjx\ + z 2 }. (17) 

In the expression for the spectral curve /z 2 = fl 2 — 2AA* understood as a polynomial 
in A, the coefficients are integrals of the motion (they are determined and given by 
(14))- 
So 

y 2 - fl 2 = -2AA* = -2(|y| 2 + (yx + xy)\ + |^| 2 A 2 ) 
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is a known polynomial and we can determine \x\ 2 , \y\ 2 and argy — aigx = arg (|). 

The previous theorem gives us a possibillity for integration of the Apel'rot case 
if we know the motion of the point v. As we will show, the motion of the point v 
is linearized on the Jacobian Jac(r 2 ) = T 2 . 

In all propositions untill the end of this section, we will assume that the Apel'rot 
conditions are satisfied. 

Theorem 3. If v\ and are coordinates of point v e T 2 , then: 

d 1 

dt vx = TS» 

Lemma 2. The functions x and y satisfy the system 

dx ( 1 1 \ iMi rz ; 

zx + i^xf) + z^y 

(18) 



dt \h I 2 ) ay/2 

dy _ ijaji + ^73) ^ iMi f I 1_ 

dt ~ I 2 X as/2 \h h 

Proof. From the Apelrot conditions (2) we can get 

/3 2 / 1 1 \ 1 1 



,{aM 2 - iM 3 ) + i^Jxl + zly 



J 1 h J I2 h 
From (19), (3), and the identity a 2 + (i 2 = 1 one obtaines 
dx _ ( 1 1 \ Mj 
dt ~ \h ~ h) ^71 v 

dy _ _ 6(071 +^73) iMi f 1 1 

dt ~ h X ay/2\h h) V 

The claim follows applying (3). 

As a consequence we get 

Lemma 3. The function v\ = — | satisfies the equation: 

dt V x) ~ x 2 \dt V dt X ) ~ x 2 



(19) 



h 

Proof of theorem 3. The divisor v is given by v\ = — | , = — il \\=-v ■ From 
lemma 2 and lemma 3, we have: 

vx ~dt\ x)- x 2 \dt v x dt - x 2 i v x ° +z ° y + 



dt 



h 



T 2 (" n !*="*) = Tfr 

(20) 
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Since T 2 is a curve of genus 1, there is only one linearly independent holomorphic 
differential 

dX _ dX 

M ~~ - 2AA* ' 

on it. 

So any linear flow on Jac(r2) is of the form = const ■ fj,. 
A consequence of (20) is that the flow of the point v is linear on Jac(T2). We 
have the next theorem: 

Theorem 4. The integration of the motion of the point v reduces to the inversion 
of the elliptical integral 



f 



d\ 1 
— t. 



,„ V^ 2 - 2AA* h 
From theorem 2 and expression (10) for x and y, it follows that, if we know the 
motion of the point v 1 in order to determine the matrix L(X) as a function of time, 
we need argai as a function of time. Let us denote by <j> x — argx, and u = tg^f. 

Theorem 5. The function u(t) satisfies the Riccati equation: 

^ = [(f(t)-9(t)]u 2 + [f(t)+g(t)}, (21) 



where 



K_ _ M _ Q\x\ 
2\x\ 



/(t) = ^i5 9(t) = 



K = Vjii^ (22) 



a \12 IiJ 

and \x\ is a known function of time. 

Proof. From 



Iim —Mi 

tgcpx 



Re x (3M X - aM 3 ' 
by differentiating and using the equations of motion (4), we get: 

d(f> x 1 = n ' ' aKl2 IlJ y^I+^g 
dt cos 2 4> x ~ (M, M) - M| 



(23) 
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From the system 

aM 1 + (3M 3 = 

(3Mx - aM 3 = \/2Rcx = V2\x\ cos<p x , 
Mi can be expressed as a function of \x\, and (f> x 

Mi = \/2(3\x\ cos<p x . 

Also, from (10) we have: 

M 2 = -V2lmx = -V2\x\ sin^. 

From (10), using (3) wc have that 2|ir| 2 = (M,M). Substituting the last three 
expressions into (23) gives 

d(j) x _ K - Q\x\ 3 cos(j) x 

By the change of variables u = tg^f the equation (24) takes the form (21). 

The classical integration procedure in the so-called Hesse coordinates [5] also 
yields the Riccati equation (see [5]). By Nekrasov (see [5], [8]) it was reduced to a 
second order linear differential equation with double-periodical coefficients. 



4. A new integrable case on so(4) x so(4) 

The equations of rotations of a heavy n-dimensional rigid body fixed at a point 
on so(n) x so(n) are given in [9]. In the moving frame, these equations are: 

^ = [M,Q] + [T,X], f = [I\ft]. (25) 

M = in + £11 <G so(n) is the kinetic momentum, £1 e so(n) is the angular velocity, 
/ is a symmetric n x n matrix, T G so(n), and X G so(n) is a given constant 
matrix. By choosing the moving frame with / diagonal, / = diag(Ii, /„), wc 
have Mij = (7, + and Jj + Ij are the principal inertia momenta. 

In [9], the Lagrange case was defined by I\ = I2 = a; I3 = ••• = /„ = fe, X i2 ^ 
0, Xij =0, i, j 7^ 1, 2, i < j. The completely symmetric case was defined there by 
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Ii = ■ ■ ■ = I n = a, where X e so(n) is an arbitrary constant matrix. It was shown 
in [9] that the equations (25) in these cases could be represented by the following 
L-A pair: 

^(A 2 C + AM + T) = [A 2 C + AM + T, XX + Q] , (26) 

where in the Lagrange case C = (a + b)X, and in the symmetric case C = aX . 

The theorem 4.1 in [9], states that there are no other cases of n-dimensional rigid 
bodies, with the equations equivalent to (26). 

However, let us consider a 4-dimensional rigid body with the conditions: 

h = h = a, I 3 = h = b, 

(27) 

AWO, X^^O, X ij= 0, z,j^{(1,2),(3,4)}, i < j. 

Note that this is neither the Lagrange nor the completely symmetrical case. The 
next theorem shows that this example satisfies all assumptions of the theorem 4.1 
from [9]. So, the system (27) is a counterexample to that theorem. 

Theorem 6. The equations (25) for a 4~ dimensional rigid body which satisfies (27) 
are equivalent to (26) with the matrix C = (a + b)X . 

Proof. We have [C, X] = from (27). Also, [C, fi] + [M, X] = 0, so (26) is equivalent 
to (25). 

One can naturaly consider (25) for n — 4 as an equation on the semidirect 
product so(4) x so(4) with the following Poisson structure on the orbits of the 
coadjoint action (see [9]): 

{/,<?}(m,^) = -fi([dif(fi,f),d 1 g(n,i')]) 

- v{[d 1 f{^,v),d 2 g{^,v)\) 

- v{[d2f{lJ.,i'),d 1 g{iJi,,v)]), 

where /, g are restrictions of the functions / and g to the orbit of the coadjoint 
action, and dif are the partial derivatives of df . From the L-A pair representation 
(26), we have that integrals of the motion for the system given by (27), are the 
coefficients in the A polynomials Tr ( A 2 C + AM + T) 2 , and (Dct(A 2 C + AM + T)) 1 / 2 . 



12 VLADIMIR DRAGOVIC, BORISLAV GAJIC 

Using the same arguments as in [9], it can be proved that the system given by (27) 
has the following four Casimir functions: 

Ji = m 34 t 12 + Mi 2 r 34 + Mi 4 r 23 + M 23 r 14 - m 24 t 13 - Mi 3 r 24 

■h = r 34 ri 2 + r 23 ri 4 — ri 3 r 24 

J 3 = r 2 2 + v\ 3 + r 2 4 + r 23 + r 24 + r 34 

Ji = Mi 2 r 12 + Afi 3 ri 3 + Mi 4 r 14 + M 23 r 23 + Af 24 r 24 + M 34 r 34 . 

Also, it follows that this system has four integrals: 
F 1 = C 12 M 12 + C 34 M 34 
F 2 = C 34 Mi 2 + Ci 2 M 34 

F 3 = M 12 M 3i + M 23 M 14 - Mi 3 M 24 + C 34 ri 2 + Ci 2 r 34 

F 4 = Mf 2 + M\ 3 + + M| 3 + M 2 2 4 + M 3 2 4 + 2C 12 T 12 + 2C 34 r 34 , 

which arc in involution. Thus, we have 

Proposition. The system (27) is completely integrable. 

Except the three cases mentioned above, there are no other systems of the type 
(25) which are equivalent to (26). This is proved in the following: 

Theorem 7. Let us suppose that X\ 2 ^ 0. The Euler-Poisson equations (25) can 
be written in the form (26) if and only if the equations (25) describe: 

a) for n — 4, the motion of the Lagrange top, the completely symmetric top or a 
rigid body which satisfies (27). 

b) for b^4, the motion of the Lagrange or completely symmetric top. 

Proof. The equations (26) are equivalent to (25), if and only if [C, X] = and 
[C, 0] + [M, X] = 0. The second relation is equivalent to 

C ij = (I i + I k )X ij , C ij = (I j +I k )X ij , k^i,j. (28) 

From C\ 2 ^ and (28), for i = l,j = 2, it follows that Ii = I 2 , L k = Ii for all 
k, I 7^ 1,2. Let us fix i, j different from 1 and 2. From (28), we have Ii = Ij ,I k = h, 
for all k, I ^ For n^4, this means that I\ = ■ ■ ■ = I n and X is an arbitrary 
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matrix, or I\ = I2 = a, I3 = ■ ■ ■ = I n and Xy = for i, j 7^ 1, 2. For n = 4, wc 
get one more case which satisfies the conditions (27). For such dj, the relation 
[C, A"] = is satisfied. 
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